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00 ' Abstract. Recent numerical simulations of radiatively inefficient accretion flows onto 

I compact objects have shown that convection is a general feature in such flows. Dis- 

sipation of rotational and gravitational energies in the accretion flows results in in- 
, ward increase of entropy and development of efficient convective motions. Convection- 

' dominated accretion flows (CDAFs) have a structure that is modified significantly in 

I comparison with the canonical advection-dominated and Bondi-like accretion flows. 

' The flows are characterized by the flattened radial density profiles, p{R) oc 

and have reduced mass accretion rates. Convection transports outward a significant 
amount of the released binding energy of the accretion fiow. We discuss basic dynam- 
ical and observational properties of CDAFs using numerical models and self-similar 
analytical solutions. 



S INTRODUCTION 

^ ' Observations of accreting black holes of different mass, from the stellar mass 

black holes to the supermassive black holes in the center of galaxies, show impres- 
^ ■ sive similarities of data, which point at identical physical processes in accreting 

plasma. Also, black hole candidates demonstrate a great variety of physical condi- 
tions in the flows and, possibly, existence of a variety of accretion regimes. Existing 
theories describe different regimes of black hole accretion flows, which can be re- 
alized under different physical conditions. If matter accretes with a low specific 
angular momentum, j <C RgC, it forms spherical flows, which are described by 
Bondi solution [1] in the case of adiabatic flows. Here Rg is the gravitational ra- 
dius of the black hole, and c is the speed of light. If matter has a large specific 
angular momentum, j ^ RgC, then accretion disks are formed. The structure of 
accretion disks crucially depends on the efficiency of radiative cooling. If matter 
radiates efficiently, i.e. the radiative cooling time is shorter than the accretion 
time, trad ^ taccr, the disks are geometrically thin, H <^ R, where H is the scale 
height of the flow, and R is the corresponding radius. The radial structure of such 
disks is described by the Shakura & Sunyaev model [2]. In the case of radiatively 



inefficient fiows, trad ^ taccr, tlie internal energy of accreting matter is close to 
the virial energy, and the formed disks are geometrically thick, H R. The thick 
disks can be formed in two limiting cases of very high, M ^ LeiIci/c'^, and very 
low, M <^ LEdd/c^, accretion rates, where LEdd is the Eddington luminosity of the 
black hole. In the former case, the optical depth of accretion flow is very large, and 
photons, which carry most of the internal energy, are trapped inside the inflowing 
matter and can not be radiated away. In the latter case, the accreting plasma is 
very diluted, so that different radiative mechanisms are inefficient to cool plasma on 
the accretion time scales. A model which describes the geometrically thick accre- 
tion disks was called the advection-dominated accretion flow (ADAF) and attracted 
a considerable attention during the last two decades [3-10]. The main feature of 
ADAFs is that most of the locally released gravitational and rotational energies 
of accretion flow is advected inward in the form of the gas internal energy, and 
the latter is flnally absorbed by the black hole. The recent interest to ADAFs was 
generated mostly by attempts using this model to explain the low-luminosity X-ray 
objects of both galactic and intergalactic nature (for recent reviews see [11,12]). 



ADVECTION-DOMINATED FLOWS 

Properties of ADAFs could be better understood by analyzing the self-similar 
solution of the height-integrated hydrodynamical equations [10]. The solution de- 
pends on two parameters, the viscosity parameter a and the adiabatic index 7, and 
satisfles the following scalings for the density p, the radial velocity vr, the angular 
velocity fl, the isothermal sound speed and H: 

p{R) oc 

vr{R) oc R-^'\ 

n{R) oc (1) 

Cs{R) oc R-^/^, 
H{R) - Cs/VLk oc R, 

where Vl^ is the Keplerian angular velocity. Because of neglection of the radiative 
energy losses of the flow, the flow structure is independent of the mass accretion 
rate M. The accretion rate determines only scales of density through the relation 
M/{At:RHvr). 

Two inconsistencies were found related to the model of ADAFs. The first is 
that ADAFs must be convectively unstable [4,5,7]. The second is connected to the 
positiveness of "Bernoulli parameter" in ADAFs [9], 

vl GM 
Be=^ + W-^>0, 



where W is the specific enthalpy. A hypothesis was proposed that the latter prob- 
lem could be solved by assuming the presence of powerful bipolar outfiows in the 

accretion flows [13,14]. However, according to the later investigations [15,16], the 
hypothesis met some difficulties. The consequences of the earlier mentioned incon- 
sistencies for ADAFs was understood after performing two- and three-dimensional 
hydrodynamic simulations of the inefficiently radiative accretion flows [17-22] . The 
simulations have revealed a new accretion regime of inefficiently radiated plasma, in 
which convection plays a dominant role determining structure and dynamics of the 
flow. Such a regime we shall refer as CDAF, the convection-dominated accretion 
flow. 

CONVECTION INSTABILITY 

Behavior of convcctive blobs in non-rotating medium depends sensitively on the 
superadiabatic gradient 

which determines the Brunt- Vaisala frequency N given by 

C 

where g^fj = —{l/p){dP/dR) is the radial effective gravity. When A^^ is positive, 
perturbations of the blobs have an oscillatory behavior with the frequency A^ and 
the medium is convectively stable. However, when N'^ is negative, perturbations 
have a runaway growth, leading to convection. Convection is present whenever 
AVc^ is positive, i.e. when the entropy gradient, given by 

dR 7 — 1 

is negative. This is the well-known Schwarzschild criterion. 

When there is rotation, a new frequency enters the problem and convection is no 
longer determined purely by the Brunt- Vaisala frequency. In this case, the effective 
frequency A^e// of convective blobs is given by 

where n is the epicyclic frequency; for O oc R~^^^, we have k = Q. Again, when 
N^ff is negative, the rotating flows are convectively unstable. 

When convection is developed, the convective blobs carry outward some amount 
of thermal energy. Under certain approximations (see details in [24]) the convective 
energy flux can be expressed in the form. 



ds 

^conv ^convPT , (2) 

where s is the specific entropy, T is the temperature, Vconv is the diffusion coefficient, 
defined by Vconv ~ {^'mI ^){~^lffY^'^ ■, ^"^^ is the characteristic mixing-length. 
The coefficient i^conv can also be expressed in the a-parametcrization form, Pconv = 
otcanvcll^K-i whcrc ttcon?; is a dimcnsionlcss parameter that describes the strength 
of convective diffusion; this parameter is similar to the usual Shakura & Sunyaev 
a. 

In the case of the self-similar ADAF solution (1), calculations lead to the con- 
clusion that -/Vg// < 0) i-6- ADAFs are convectively unstable and must have the 
outward energy fluxes defined by equation (2). 



NUMERICAL RESULTS 

To describe the results of numerical simulations of radiativcly inefficient accre- 
tion flows, we shall mainly follow [21,22]. In these works the accretion flows were 
simulated by solving the nonrelativistic, time-dependent Navier-Stokes equations 
with thermal conduction, 

|+pVv = 0, 

p'^ = -VP + pv* + vn, (3) 

de 

p— = -PVv- Vq + g, 

where e is the speciflc internal energy, $ is the gravitational potential of the black 
hole, n is the viscous stress tensor with all components included, q is the heat flux 
density due to thermal conduction and Q is the dissipation function. No radiative 
cooling was assumed and the ideal gas equation of state, P = (7 — l)pe, was 
adopted. Only the shear viscosity was considered, 

u = acll^K. (4) 

where a is a constant, < a < 1. 

In the simulations, it was assumed that mass is steadily injected into the calcula- 
tion domain from an equatorial torus near the outer boundary. Matter is injected 
with almost Keplerian angular momentum. Owing to viscous spread, a part of the 
injected matter moves inward and forms an accretion flow. The computations were 
started from an initial state in which there is a small amount of mass in the do- 
main. After a time comparable to the viscous timescale, the accretion flow achieves 
a quasi-stationary behavior. 

Results obtained by using two-dimensional axisymmetric simulations are sum- 
marized in Figure 1. Four types of accretion flows can be distinguished from the 
models with no thermal conduction (see Figure 1, left panel). 



1. Convective flows. For a very smaU viscosity, a < 0.03, accretion flows are 
convectively unstable. Axially symmetric convection transports the angular mo- 
mentum inward rather than outward. Convection governs the flow structure, mak- 
ing a flattened density proflle, p{R) oc R^^^"^, with respect to the one for ADAFs. 
Convection transports a signiflcant amount (up to ~ 1%) of the dissipated binding 
energy outward. No powerful outflows are present. 

2. Large-scale circulations. For a larger, but still small viscosity, a ~ 0.1, 
accretion flows could be both stable or unstable convectively, depending on a and 
7. The flow pattern consists of the large-scale (~ R) meridional circulations. No 
powerful unbound outflows are presented. In some respect this type of flow is the 
limiting case of the convective flows in which the small-scale motions are suppressed 
by larger viscosity. 

3. Pure inflows. With an increasing viscosity, a ~ 0.3, the convective instability 
dies off. Some models (with 7 ~ 3/2) are characterized by a pure inflow pattern, 
and agree in many respects with the self-similar ADAFs. No outflows are present. 

4. Bipolar outflows. For a large viscosity, a ~ 1, accretion flows differ consid- 
erably from the simple self-similar models. Powerful unbound bipolar outflows are 
present. 
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FIGURE 1. Properties of two-dimensional axisymmctric models of accretion flows. Models with 
no thermal conduction arc shown in the left panel. Models with thermal conduction (Prandtl num- 
ber Pr = 1) are shown in the right panel. Each circle represents a model in the (a, 7) parameter 
space. The empty circles correspond to laminar flows, the crossed circles represent unstable mod- 
els with large-scale (~ R) meridional circulations and solid circles indicate convective models. 
Two outward directed arrows correspond to bipolar outflows, whereas one arrow corresponds to 
a unipolar outflow. Two inward directed arrows correspond to a pure inflow. The dashed line on 
the left-hand panel approximately separates regions of convectively stable/unstable flows. 



Figure 2 illustrates the density distribution and flow pattern in the low viscos- 
ity models. The models show complicated time-dependent behavior with numer- 




tSi O Ci' 

in o; 



FIGURE 2. Snapshots of density distribution (upper panel) and streamlines (lower panel) from 
two-dimensional model of CDAF with a = 0.01 and 7 = 5/3, in the meridional cross-section. 
The black hole is located at the origin. The flow pattern is highly time dependent and consists 
of numerous vortices of different spatial scales. The variability of flow pattern is accompanied by 
density variations and results in variability of the accretion rate. 



ous vortices and circulations, and with density fluctuations. However, the time- 
averaged flow patterns are smooth and do not demonstrate smaU-scale features. 

Figure 3 shows typical behavior of radial proflles of variety of time- and angle- 
averaged variables in convective models. Except near the inner and outer bound- 
aries (i? < lORg and R > 10^ Rg), the proflles of p, Cs, v and Q show prominent 
power-law behaviors. 

Two-dimensional models have demonstrated that the (r,(/7)-component of the 
volume averaged Reynolds stress tensor takes negative values in the case of con- 
vective flows. This means that axisymmetric convection turbulence transports the 




FIGURE 3. Selected properties of two-dimensional axisymmetric model with a = 0.01 and 
7 = 5/3 [23]. All quantities have been averaged over polar angle 9 and time. Except near the 
boundaries {R < lORg and R > lO^Rg), the profiles of p, (?g/v%, and Q./Q,k show the power-law 
behaviors predicted by the self-similar convective envelope solution. In the upper left panel the 
dashed line corresponds to the analytical scaling p oc R~^/'^. Similarly, in the lower left panel the 
dashed line corresponds to the expected scaling v oc 

i?-3/2 for CDAFs. 



angular momentum inward. Indeed, if there are no azimuthal gradients of pressure, 
the turbulence attempts to erase the angular momentum gradient, which for the 
Keplerian-like angular velocity, oc i?"^/^, means that the angular momentum 
transports inward. However, this is a property of axisymmetric flows. Does con- 
vection in real three-dimensional flows move angular momentum inward? Do such 
flows have significant azimuthal perturbations which destroy axisymmetry? The 
answer to this question has been found with the help of fully three-dimensional 
simulations, without the limiting assumption of axisymmetry [21]. The simula- 
tions have demonstrated a good qualitative and quantitative coincidence of the 
results of two- and three-dimension modeling. It was demonstrated that in the 
differentially rotating flows with nearly Keplerian angular velocities, all azimuthal 
gradients are efficiently washed out, leading to almost axisymmetrical structure of 
three-dimensional flows. Figure 4 illustrates axisymmetric structure of the low- 
viscosity accretion flow obtained in three-dimensional simulations. 

All low-viscosity models have negative volume-averaged Be in all ranges of the 




FIGURE 4. Snapshot of surfaces of a constant density in three-dimensional model of CDAF 
with a — 0.01 and 7 = 5/3. Only a quarter of the innermost region of the full domain is shown. 
Accretion matter rotates around 2-axis. The black hole is located at the origin. All axes are 
labeled in the units of Rg. One can clearly see that density perturbations, associated with motion 
of convective blobs in the accretion flow, have axially symmetric form with respect to the axis of 
rotation. 



radii. Only temporal convective blobs and narrow regions at the disk surfaces show 
positive Be. It seems that the outward directed convection energy transport pro- 
vides an additional cooling mechanism of the flows that always results in averaged 
Be < 0, opposite to positive Be obtained in ADAFs. 

In several axisymmetric simulations, effects of turbulent thermal conduction have 
been studied. It was found that the conduction has an important influence on the 
flow structure, but it docs not introduce a new type of flow in addition to those 
already discussed (sec Figure 1, right panel). The conduction leads to a suppression 
of small-scale convection in the low-viscosity models. In the moderate-viscosity 
models the thermal conduction acts as a cooling agent in the outflows, reducing or 
even suppressing them. 

Obtained properties of CDAF can have important implications for the spec- 
tra and luminosities of accreting black holes [20,22,23]. Indeed, since p oc 
and T oc R~^, the bremsstrahlung cooling rate per unit volume varies as Qhr oc 
p2j'-i/2 ^ i?-3/2 After integration of Qbr over bulk of the flow one obtains the 
bremsstrahlung luminosity L^r oc R. This means that CDAFs mostly radiate on the 
outside, whereas most of the radiative energy of ADAFs comes from the innermost 
region [10]. 



SELF-SIMILAR CONVECTIVE ENVELOPE 

Numerical simulations of convective accretion flows have been understood quali- 
tatively after construction of analytical self-similar solutions of such types of flows 
[24,25]. In the presence of convection, the height-integrated angular momentum 
and energy equations can be written as follows, 

J adv ~l~ J vise ~l~ J conv 0; (5) 



P-^-^^^ ~l~ (j^J^^^ Pconv^ ^{^Jvise ~l~ Jconv^i (6) 

where Jadvi J vise a-^d J conv are the advective, viscous and convective angular mo- 
mentum fluxes respectively. Solutions of equations (5) and (6), together with the 
equation of motion in the radial direction, crucially depend on the direction of 
convective angular momentum transport. There are two extreme possibilities. The 
first one is that convection behaves like normal viscosity, 

Jconv — —^eonvPR?{d^/dR)- (7) 

An alternative possibility is that the convective flux is directed down the speciflc 
angular momentum gradient. 



Jeonv — —^ccmvPR[d{P'R?)/dR]. 



(8) 



For Q oc convective fluxes (7) and (8) correspond to outward and inward 

transport of angular momentum respectively. Numerical simulations have demon- 
strated that Jconv < in CDAFs. 

Choosing the convective flux (8) one can construct a nonaccreting solution with 
vr = 0. We refer to it as a "convective envelope". For this solution, from equation 
(5) one obtains Jadv — and Jconv — —Jvisc- The latter means that in the convective 
envelopes the net angular momentum flux vanishes: the inward convective flux is 
exactly balanced by the outward viscous flux. Because of absence of advection 
(due to Vr = 0) and zero local dissipation rate (due to zero net angular momentum 
flux), the energy equation (6) is satisfied trivially with Fconv{R) ^ R"^- Finally, 
the convective envelope solution has the following scalings: 

p{R) oc 

vr{R) = 0, 

n{R) oc R-^^\ 

Cs(R) oc i?-^/^ 
H{R) oc R, 

The outward energy flux Fconv corresponds to luminosity Lconv — 4:Tr RHFconv, 
which is independent of radius. Source of this energy flux is formally located at 
R = 0, but the nature of the source is not specifled in the self-similar solution. In 
more realistic CDAFs, the mass accretion rate M is small, but it is not exactly 
zero. The fraction of the binding energy of the accreted mass released in the 
innermost region is the source of the energy required to support convection. In the 
case of CDAFs the convective luminosity can be expressed in the following form, 
Lconv = sMc^, where the parameter e has been estimated in numerical simulations, 
e ~ 0.01. A non-zero M leads to a finite vr = M/{At:RHp) oc R'^l^. The latter 
scaling and the scalings of other quantities in (9) agree quite well with the numerical 
results (see Figure 3). 



(9) 



RECENT MHD SIMULATIONS 

There is good reason to believe that "viscosity" in differentially rotating accretion 
fiows is produced by magnetic stress generated by the magneto-rotational instability 
(MRI) [26]. MHD simulations of accretion fiows under conditions, in which CDAFs 
could be formed, have been recently performed ([27-32] and Matsumoto et al. in 
this Proceedings). A general conclusion of these studies is that MRI, and probably 
other instabilities of magnetized medium, leads to development of turbulence in 
accretion fiows. The turbulence re-distributes angular momentum in the fiows 
with an effective Oe// ~ 0.01 — 0.1. The simulations show a tendency for density 
profiles in accretion fiows to be fiatter, than those in ADAFs. 



In general, these results of MHD simulations do not contradict the results of 
hydrodynamical simulations of convective flows, which were discussed in previous 

sections. The MHD simulations have demonstrated the range of aeff, in which 
hydrodynamical models are convectively unstable. The flattened density profiles 
in MHD models are in good agreement with those found in convective hydrody- 
namical models. However, a direct comparison of results of both approaches meets 
difficulties. The main problem is a low spatial resolution, and as a consequence, a 
small radial range of flow patterns studied in the discussed MHD simulations. In 
this case of small radial ranges the effects of boundaries become signiflcant enough 
to influence the flow structure. Thus, more powerful computers and more devel- 
oped numerical codes are required to obtain a comparable spatial resolution to that 
used in hydrodynamical simulations. 

Another problem which can introduce difficulties for the comparison is that some 
MHD simulations have uncontrolled losses of energy due to numerical reconnection 
of magnetic lines. These energy losses can artificially suppress convection in the 
considered flows and change the flow structures, independently of the used spatial 
resolution. The problem could be solved by using an "artiflcial" resistivity in MHD 
schemes (see [32] for detailed discussion). 

CONCLUSIONS 

• Radiatively inefficient hydrodynamical black hole accretion flows with small 
viscosity (a ^ 0.1) are always dominated by convection. Thus for low viscosity 
flows the concept of ADAF is unphysical. 

• At present, CDAFs provide the best theoretical model for understanding of 
radiatively inefficient, low viscosity accretion flows. 

• CDAFs are hot, their thermal energy is close to the virial energy. They have 
relatively reduced mass accretion rate. The density distribution in CDAFs is 
flattened, comparing to the one in ADAFs and Bondi-like flows (p oc 

in the former case vs. p oc in the latter case). Convection transports 

outward a signiflcant amount (up to ~ 1%) of dissipated binding energy of 
accretion flows. No powerful outflows are present in CDAFs. 

• CDAF is a very simple model which is now undergoing a rapid development. 
A more detailed account of plasma, magnetic, and radiative processes is likely 
to change this model. 
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